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Missing completely of the CMB quadrupole in WMAP data 

Hao LIU^ and Ti-Pei LI^'^ 

cn ■ ABSTRACT 

o 

, In cosmic microwave background (CMB) experiments, foreground-cleaned 

c^ . temperature maps are still contaminated by the residual dipole due to uncer- 

tainties of the Doppler dipole direction and microwave radiometer sidelobe. To 
obtain reliable CMB maps, such contamination has to be carefully removed from 
observed data. We have previously built a software package for map-making, 
O '■ residual dipole-contamination removal, and power spectrum estimation from the 

w i Wilkinson Microwave Anisotropy Probe (WMAP) raw data. This software has 

"-^ ! now been significantly improved. With the improved software, we obtain a neg- 

I , ative result of CMB quadrupole detection with the WMAP raw data, which is 

U ; -3.2 ± 3.5 fiK^ from the seven-year WMAP (WMAP7) data. This result is 

52 ! evidently incompatible with ~ 1000 /xK^ expected from the standard cosmolog- 

ical model ACDM. The completely missing of CMB quadrupole poses a serious 
CN ! challenge to the standard model and sets a strong constraint on possible models 

Q i of cosmology. Due to the importance of this result for understanding the origin 

CN I and early evolution of our universe, the software codes we used are opened for 

i/-^ [ public checking. 

rn 

i:^ ! Subject headings: cosmology — cosmic microwave background — data analysis 



/\ . 1. Introduction 

The CMB quadrupole is the largest observable structure in our universe. It corresponds 
to the temperature fluctuation power at the lowest detectable multiple moment / = 2, or ~ 
90 ° angular scale, and reflects the universe circumstance at very early epoch. It is a difficult 
task to truly evaluate the quadrupole moment from an observed microwave map. We find 
that, although the foreground contamination problem has been effectively treated in CMB 
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observations, another important systematic error, the scan-induced dipole contamination 
accumulated on large angular scales, s till exists and produces a false aligned quadrupole 
for all sweep missions ( iLiu fc Lill2011al Jbl). which, unfortunately, has not been adequately 
realized and removed from the official results released by the CMB mission groups. 



For a sky pixel with the unit direction vector n, the Doppler dipole signal is 

To 



T. 



-V ■ n. 



1^ 



where Tq = 2.725 K is the monopole temperature, v the joint velocity of antenna relative 
to the CMB, c the speed of light. The amplitude of Doppler dipole moment is greater 
than 3000 /iK, much stronger than the CMB anisotropy (~ 50 /iK). The dipole signal T^ 
must be calculated with Eq.[T] for each observation and removed from the measured tem- 
perature. However, the uncertainties of v and n in Eq.[T] make it impossible to estimate 
Tj acc urate enough. Firstly , the direction error ~ 10' in the vector v, i.e. in the CMB 
dipole (JBennett et al.ll2003al ). can cause the calculated dipole signal to deviate by as much as 
10 — 20/xK, which cannot be ignored when compared to the relatively weak CMB anisotropy. 
Secondly, in calculating T^ with Eq.dl the effect of overall sidelobe response uncertainty of 
radio telescope can equivalently in troduce an erro r in the direction vector n, which is esti- 
mated to be ~ 10'— 20' for WMAP fJLiu &: Lill2011al ). Furthermore, an undoc umented 25.6 ms 
timing offset error in WMAP time-ordered-data (TOD) has been revealed (JLiu. Xiong &: Li 
20 111 : lRoukemall2010l ). which also can deviate calculated dipole signals as an antenna point- 
ing direction error of ~ 7' does . Different errors in dipole calculation with Eq.[T] can be 
synthesized to an overall-equivalent direction error An with a magnitude of about 10'. Since 
An includes both real and equivalent pointing error, it can not be eliminated by improving 
only the real antenna pointing accuracy. After a full sky scan, these An induced dipole de- 
viations should be accumulated in the resulting CMB map, producing artificial anisotropics 
at the large angular scales (not only the dipol e) with a pattern closely correlated to the scan 
pattern and an artificial aligned quadrupole (JLiu fc Lill2011bl ). 



The WMAP scan scheme makes its sky coverage inhomogeneous - the number of obser- 
vations being greatest at the ecliptic poles and the plane being most sparsely observed. 
As above expected, it is really existed in released WMAP temperature maps a signifi- 
cant large-scale non-Gau ssian modulation feature being closely correlated to the scan pat- 
tern (ISpergel et al.l 120061 ) . and the CMB quadrupole is oddly aligned along the scan pattern 



(see the first plot of Fig. [TD ca lled the "axis of evil" (Ide Oliveira-Costall2004l : ISchwaztz et al. 



2004 : iLand fc Magueijd l2005l ) . With an assumed overall-equivalent direction error An of 
7', we simulated the first year WMAP scan survey, calculated the dipole deviation for each 
observation, and found that the accumulated temperature distortion, which is generated 
using only spacecraft attitude information and the scan scheme without any CMB informa- 
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tion a t all, is well consistent with the CMB quadrupole component released by the WMAP 



team (JLiu. Xiong &: Lill201l[ ). Therefore, without removing the accumulated dipole contam- 



ination, released CMB anisotropy at the largest scales is not reliable. 

To give the exact overall-equivalent direction error An by estimating all effects caus- 
ing the dipole contamination is impossible, because the sidelobe uncertainty can not be 
ehminated completely. Independent of the reasons that caus e An, we have pr oposed a tem- 



plate fitting procedure to remove the dipole contamination (JLiu fc Lill2011al ). The overall- 
equivalent direction error can be expressed as An = Ax -|- Ay -|- Az in the spacecraft 
coordinate system {X,Y,Z), where the X-axis is parallel to plane of radiators, the Z-axis 
is the anti-sun direction of the spin axis, and the F-axis is perpendicular to both. For an 
assumed overall direction error An = Ax = 1', we can follow the WMAP observational 
scan history to calculate dipole deviations by ST^ = (To/c)v ■ An for each observation and 
accumulate them into a template map AT^ . The other two template maps ATy and AT^ can 
be similarly derived by assuming An of 1' being along the Y- and Z-axis, respectively. With 
template cleaning, from an observed and monopole subtracted map T' the clean temperature 
map T can be derived by 

T = T'- {a^AT, + ayAty + a,At,) , (2) 

where the coefficients Qx-, ay and a^ are determined by minimizing the variance of T, 

T^ = min . (3) 

The template maps ATj., ATy and AT^ have to be generated from the WMAP original 
raw data with map-making software. To remove the accumulated dipole contamination with 
Eq. [2] and [3] to produce cleaned CMB maps, we must redo the whole map-making process 
from WMAP raw TOElJ. Because the WMAP map-making software has not been opened for 
public use, we built a self-consistent software package of map-making and power spectrum 
estimation. With our previous software, we modeled and removed the pseudo-dipole signal 
from the released WMAP 7 maps of the Q-, V- and W-bands and then from the clean maps 
obtained the average residual quadr upole power of ^ 17.1 /xK^, only ~ 14% of what has 



been released by the WMAP team ( iLiu fc Lil l2011al ). Since then our software has been 



much improved, and in this work we use the improved software and obtain the amplitude 
of CMB quadrupole as —3.2 ± 3.5 /xK^ from the WMAP7 data, which indicates that the 
CMB quadrupole predicted by the standard cosmology model might be completely missing. 



^Because the deviation templates iSTx^y^z must be calculated with the same pipeline for the temperature 
map-making. 
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2. Analysis software 



Independently of the WMAP team, we developed a self-consistent softwar e package 



of map -making: an d power spectrura estimation, which passed a variety of tests (JLiu fc Li 



2009af ) and used in lLiu fc Lil fl2011af ). 



2.1. Improvements 



The software we used in iLiu &: Lil (J2011al ) has many deficiencies as briefiy listed below: 



» The antenna pointing error is a preset absolute value. 

» It does not distinguish the difference between real and equivalent pointing errors. 

» The temperatures inside the processing mask are not calculated. 

» The foreground removal and power spectrum calculation have not been included as a 
whole package. 

» It is not easy to use. 

Now the software has been much improved: 

» The antenna pointing error is determined by template fitting, no preset value any more. 

» The templates are calculated for each single year and single band separately. 

» We consider mainly the equivalent pointing error now. 

» The correct loss- imbalance factors are given automatically according to the TOD ver- 
sion (WMAP3, 5, 7). 

» The produced maps are now full-sky. 

» The software now starts from the TOD and ends at the final CMB cross-power spec- 
trum. 

» Now all intermediate steps are done automatically by typing just a short command. 
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It is important to note that, all improvements to the software are "fair" for quadrupole 
detection. This means that, in principle, they do not prefer a higher or lower quadrupole by 
themselves. This will be further discussed in §|11 

The current software contains the following modules, which will be discussed one by 
one: 

1. Unpack the TOD 

2. Make the CMB temperature maps 

3. Remove the foreground emission 

4. Remove the large scale deviations 

5. Compute the output CMB power spectrum 

2.2. Modules in this version 

The first module is to unpack the TOD. In this module, we compute the A- and B-side 
antenna pointing vectors n^ and n^ , the differential Doppler signal 

rf = -v-(n,-nj (4) 

c 

with To being the CMB monopole, c the light speed, v the velocity of the spacecraft relative 
to the CMB rest frame, and the differential dipole deviations 

d*^=^^r■ An (5) 

with An^, An , An^ being the e quivalent pointing errors along the X, Y, Z axis of the 



spacecraft coordinate (see Eq. 4 of iLiu &: Lil (l2011al )). Since they will eventually be deter- 



mined by fitting, the absolute amplitudes of An^, An^^, An^ are not important here. Unlike 
the previous version, in this module we do not apply any data selection. All data will be 
passed to the next module. 

The second module is to make CMB temperature maps. Here we exclude some of 
the data by TOD flags and planet po si tion, and make the output CMB temperature maps 



according to Eq. 19 of iHinshaw et al.l (120031 ) with 80 rounds of iterations. The default is 
to make single year maps, because the CMB cross power spectrum should be derived from 
single year maps. The main difference to the previous version is: The data selection is done 
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here, so we do not have to re-unpack the TOD for a different data selection criterion. This 
saves a lot of time because the unpacking is the most time-consuming part of all. In this 
step, we also make three temperature deviation templates ATj., ATy and AT^ from d*, d* 
and (i*, with exactly the same condition to the corresponding single year CMB temperature 
map. The masked region is first ignored, and then, after other regions have been generated, 
they are used to re-calculate the masked region to give a full-sky output. Due to the planet 
exclusion, some places might still be empty; however, they are no more than 0.1% of the 
full sky, so we just fill them with zero to obtain a full sky output, and set the number of 
observations of those pixels to zero so that they can be easily identified. 

After ma k ing ma ps, we then smooth them and remove the foreground according to 



Bennett et al.l ( l2003bl ). In previous work, we simply use the WMAP coefficients for fore- 
ground subtraction; however, this time the coefficients are self-consistently determined by 
data, which is apparently more reliable. 

The next step is to remove AT^, ATy and ATz from the corresponding single year map 
T' by template fitting with Eqs.[2]and|n]to get the foreground cleaned and dipole deviation 
fixed CMB temperature map T. In previous version, considering the correlation between 
ATj., ATy, ATz, we use only ATj. in fitting. This time, we first remove the components 
that are proportional to AT^ from ATy and ATz, and then remove the component that is 
proportional to the remaining ATy from the remaining ATz to get three linearly independent 
templates, which are removed from T' by fitting. 

In the last step, we calculate all CM B power spectra from the cross of all single year 



CMB temperature maps using PolSpice ( ISzapudi et al.l I2OOII : IChon et al.ll2004J ). and then 



bin the average of all cross power spectra according to the WMAP binning scheme. Both 
binned and unbinned power spectra are provided as the final outputs. 



2.3. Power spectra 

Using our software we produce CMB temperature anisotropy maps T{n) from WMAP 
TOD data, and decompose T(n) in spherical harmonics 

I 

^(j^) = 5Z 5Z (^imyimin) (6) 

/>0 m=-l 
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Fig. 1. — From left to right, the first plot is the released CMB quadrupole component derived 
from WMAP5 V and W bands. The other three plots are examples of the X-, Y-, and Z- 
templates (WMAP7, year-1, Wl-band) used in our new code to remove deviations due to 
the overall-equivalent direction error. The new templates are significantly improved: They 
are whole-sky maps, and each single-year, single-band map has its own X- Y- Z- templates. 
All plots are in Galactic coordinates. 

where n is a unit direction vector. For different multiple moment / we calculate the angular 
power /(/ + l)Q/27r where 

Ci = ^^ y \a, \\ (7) 

m=—l 



3. Analysis result 

For each single year WMAP observation and each band, we calculate the accumulated 
dipole contamination maps ATj;, ATy and AT^ as removal templates with an assumed equiv- 
alent pointing error An^ = An^^ = An^ = 1'. As an example, the later three plots of Fig. [1] 
show the obtained template maps AT^, ATy and AT^ from the WMAP7 year-1 Wl-band 
data. 

Now we can remove the accumulated dipole deviation with our map-making software. 
We produce single-year CMB temperature maps from the WMAP3, WMAP5 and WMAP7 
data separately, and clean for the foreground as done by WMAP, and remove the dipole 
contamination by template fitting by Eqs.|2] and [31 After that, cross-power spectra are 
calculated from each pair of the clean maps, and the average of all cross-power spectra is 
used as the final CMB power spectrum. 

From Fig. [21 we notice that our resulting spectrum remains highly consistent with the 
WMAP release at high-/, either before or after the dipole contamination removal. This 
means that our software is fully consistent to WMAP and the deviation removal which 
affects mainly the low-/. 

Unlike the high-/, the quadrupole and octopole results are significantly different to the 
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Fig. 2. — High / CMB angular power spectra. Cross: the WMAP7 official release. Circle: the 
spectrum derived from WMAP7 data using our software without removing dipole deviation. 
Filled circle: the spectrum derived from WMAP7 data by our software after removing dipole 
deviation. 

WMAP release. The amplitudes of CMB quadrupole from WMAP3, WMAP5 and WMAP7 
data drop from the WMAP release of greater than 200 /iK^ to between —4 and —0.1 yuK^, as 
shown in Fig. [3] and Table [H The amplitude error mainly consists of two components: map- 
to-map variance (due to noise, residual foreground, etc) and incomplete sky coverage (due 
to the foreground mask). The map-to- map variance can be directly derived from the cross- 
power spectrao, and the incomplete sky coverage effect can be estimated by simulations. 
In this way, the standard error of the cleaned CMB quadrupole of WMAP7 is estimated 
to be 3.5 /iK^, very close to the standard error give by the WMAP team here (~ 5/iK^). 
Therefore, the result of CMB quadrupole detection from the seven year WMAP observation 
can be reported as Q = —3.2 ± 3.5 yuK^. 



4. Statistical significance 



4.1. The degree of freedom problem 

In stead of only one template AT^ used in iLiu fc Lil (120 11 al ibi) to remove pseudo-dipole 
signal, now we use all three templates AT^, ATy and AT^. Since the quadrupole (/ = 2) has 



^Each independent single-year, single-band temperature map will be used repeatedly in calculating the 
cross power spectra; however, this does not significantly affect the variance calculation, because each single- 
map has the same contribution in the final cross power spectra. 
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3. — Low / CMB power spectra. Filled circle: our result from WMAP7 data after the 



dipole contamination removal. Cross: the WMAP7 official release. The unit is yuK^. All the 
error bars in this figure are small than the marker symbols. 

only five degrees of freedom (see Eq. [6] and [7] where m from -2 to 2), theoretically speaking, 
if we have five linearly independent quadrupole templates, it will be possible to remove the 
quadrupole completely from any map by using Eq. O Although we have only three templates: 
ATx, ATy and ATj, one may still worry about the degree of freedom problem. 

However, it must be noticed that, the precondition for degree of freedom problem is that 
the templates must be quadrupole templates, or at least dominated by the quadrupole. This 
is true for the AT^ template, but false for ATy and ATz templates. As shown in Fig. HI both 
ATy and AT^ are dominated by non-quadrupole components. Therefore, when we add these 
two templates into fitting, the degree of freedom will not change. In other words, by adding 
these two templates into fitting, the probability distribution of the resulting quadrupole 
power will not change significantly. 



Table 1: CMB quadrupole and octopole (/xK^) 



Data 


Quadrupole 


Octopole 




WMAP release This work 


WMAP release This work 


WMAP3 


211 -0.1 


1041 430 


WMAP5 


213 -4.0 


1038 415 


WMAP7 


201 -3.2 


1050 411 



-10- 



0.010 

o.ooa - 

0.006 
0.004 



200 



lis 



400 500 




800 1000 1200 




1000 1200 



Fig. 4. — The power spectra of ATy (left) and AT^ (right) templates of year-1, Wl. The 
dipoles and monopoles have been removed in advance, and the correlated components be- 
tween AT^, ATy and AT^ have also been canceled as illustrated in § 12.21 The maximum 
power for ATy occurs at / = 8, and the maximum power for ATj occurs at / = 931. 



This can be tested by simulation: We generate 3 x 10^ simulated CMB temperature 
maps according to the best-fit WMAP CMB power spectrum, with the CMB quadrupole 
Cil{l + l)/27r (/ = 2) set to 100 /iK^, and apply the deviation template fitting with the same 
process applied to the real data, then compute the ratio R = C2/C2, where Cg is obtained 
before fitting and C2 after fitting. The distribution of R for only one template ATj. used 
is shown by the dotted line in Fig. |5l and that for three templates by the solid line. It is 
apparent that, by adding the y- and 2;-templates into fitting, the probability density function 
does not change significantly. The average of R is 0.761 for three templates condition, and 
0.813 for x-only condition, very close to each other. We can also see that for both case, many 
simulated map can give R > 1, which means the fitting even has a good chance to increase 
the quadrupole. 

We should also notice that, with either methods {xyz or x-only), the average of R will 
be systematically lower than 1. This is inevitable due to chanced correlation between the 
templates and true CMB - when we try to remove the foreground with emission templates, 
the same thing also happens. We never deny such a deficiency. However, as we can see 
from simulation, the average of R is close to 1. This tells us that, before finding a better 
method to calculate the exact accumulated dipolar error, it is still necessary to clean the 
contamination using Eqs. [2] and El and keep aware of the possible deficiency, like what we 
have done in cleaning the foreground. 
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Fig. 5. — The probability density functions oi R = C2/C2, where Cg is derived before 
template fitting, and C2 after fitting. Solid line: all x-, y- and ^-templates are used in 
fitting. Dotted line: only the x-template is used in fitting. 

4.2. What does the negative quadrupole mean? 

Another problem is the CMB quadrupole value after removing dipole contamination be- 
ing negative. However, from Eq.[7]it must be greater than or equal to zero. Incomplete sky 
coverage due to the galactic mask causes uncertainty to the final power spectrum. Therefore, 
if the true CMB quadrupole is zero or very low, we have equal chance to get positive or nega- 
tive estimations, otherwise the estimation will be biased. In other words, with the existence 
of the Galactic mask and zero/nearly zero quadrupole, negative quadrupole estimation is 
just a natural result of statistical uncertainty. 

To correctly understand the statistical meaning of a negative C2, the best way is to use 
it to find the "most likely" true quadrupole power ((72''"'^). This means: Given the CMB 
quadrupole value estimated from the experimental data set (C2, can be negative), what's 
the most likely true CMB quadrupole (6*2*^"^, must be positive or zero)? This can be done 
by using the Bayes's theorem: 



P(C*™"|C2) 



P(C2 



Where P{C2^'^'^\C2) is what we want: the conditional probability of the true quadrupole 
power given C2, and P{C2\Cl^^'^) is the conditional probability of getting C2 given a true 
CMB quadrupole, which can be obtained by simulation. 

In general case, it is quite difficult to calculate Eq. [HI however, here the problem can 
be greatly simplified: For a given data set, C2 is a constant, so do P{C2)- P{C2\C2^^) is a 
monotone decreasing function when the constant C2 < 0. As for P{C^^'^), in most cases, 
when there is no a priory limitation upon Cg*""^, it is a convention to assume PiC^'^'^) = 1. 
Now we can say that, when C2 < 0, Eq. |8]will be maximized at Cg™^ = 0. Stronger negative 
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value for C2 will only make P{Cl^'^'^\C2) decrease more quickly along the positive X axis. 
Therefore, although we can not give the exact value of P(C2''"^|C2), we can conclude that 
its most likely value of Cg*""^ is exactly zero. 



4.3. On the MCL method 

The maximum likelihood-C; (MCL) method is preferred in estimating bot h the low-/ 
CMB power spectrum an d the model-data agreements in the WMAP results (lEfstathiqu 



2004J : iHinshaw et all 120071 : iDunklev et al.l 120091 : lEfstathiou et al.l 120101 : iLarson et al.ll2011f ). 
The measured CMB quadrupole of ~ 120 fiK"^ from the WMAP7 data has been signifi cantly 
increased up to ~ 200 /iK^ after using MCL by the WMAP team jBennett et al.l I20 111 ). The 



core of all forms of MCL is the covariance matrix C =< xx^ > and its inverse matrix 
C^^. C can also be written in spherical harmonic space as Cw = 6{l,l')Ci (ignoring the 
noise, Ci is the "true" CMB power spectrum). In general cases, this is fine; however, now 
we have seen the possibility that the true CMB quadrupole might be zero, which can make 
C^^ non-existent. Another thing also need to be noticed: The MCL method can never 
give a negative result for the low-/ power spectrum, because negative result has no likelihood 
meaning. However, as we discussed in § 14.21 when there is a mask and the "true" quadrupole 
is suspected to be zero, negative estimators is necessary to ensure that the expectation of 
estimated power spectrum is unbiased (all positive estimators will only give a > average 
as expectation). Therefore the MCL method needs to be further studied. 



4.4. Significance 

To estimate the significance of the negative result of CMB quadrupole detection, we 
generate simulate d ACDM expected C MB sky maps with the synfast routine in the HEALPix 
software package (JGorski et al.l 120051 ). For each simulated map, we calculate the magnitude 
of quadrupole Q' without dipole contamination cleaning and Q after the template cleaning 
with Eq.[2] and Eq.[3l From 3 x 10^ simulated maps with Q' = 100 /iK^, we find that only 
4.1 X 10"^ of them have Q < 0.1 /iK^, and only 1.7 x 10"^ of them have Q < -3.2 /iK^. 
Therefore, the chance that the negative result of CMB quadrupole detection with WMAP 
comes from the real CMB quadrupole being occasionally identical with the accumulated 
dipole deviation is very low. 
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Discussion 



According to ACDM, the lambda dark matter model with an initial inflation, the primor- 
dial density fluctuations have a nearly scale-invariant spectrum. The CMB quadrupole rep- 
resents spatial temperature fluctuation on large angular scales of ~ 90°. From the flrst-year 
spectrum of CMB anisotropies the WMAP team derived an CM B quadrupole of ~ 2 00 /xK^, 
which is unexpectedly low than ~ 1000 /iK^ predicted by ACDM (ISpergel et al.ll2003l ). Since 
the flrst release of WMAP results, quite a few models have been proposed to explain the 
strange loss of power by attributing new physics at work in the early universe. As a notice- 
able example, a small unive rse with a specific ri gid topology was suggested to account the 
low quadrupole amplitude (JLuminet et al.ll2003l ). Whereas now after removing the dipole 
contamination the quadrupole is further missing completely in WMAP data, which is hard 
to be explained along by modifying the geometry or inflationary model of early universe and 
poses a serious challenge to current cosmological models. 

The measurement of very weak anisotropies of CMB temperature is an extremely com- 
plicated and hard task. As the flrst mission with unprecedented accuracy able to make 
the study of cosmology with high precision, it is needed time and common effort by cos- 
mology society to reveal unrealized systematic and statistical errors in the experiment step 
by step. Thanks to the WMAP team for making WMAP data publicly available, from re- 
analyzing released WMAP raw data we realize the previously ignored problem of artiflcial 
CMB anisotropies caused by coupling of dipole deviation and obs ervational scan a nd propose 
the method to remove it from observed data of CMB missions (JLiu fc Lill2011al Jbl). Due to 
the importance of the missing of CMB quadrupole for understanding the origin and early 
evolution of our universe, the software codes we used are opened for public checking and 



usm 

Besides the dipole contamination, we have discovered other kinds of systematic error in 
released CMB maps, e.g., scan-rings of hot sources being signiflcantly cooled by ~ lO^uK, 
and an up to ~ 20 /xK syste matic deviation of average temperature in Galactic latitude 
distributions (JLiu fc Lill2009al ). te mperature disto rtions up to ~ 20 /xK caused by the trans- 
mission imbalance of radiometers (JLi et al.ll2009l ). These remarkable observational effects 
should distort powers at larger multiple moments and have to be removed from the CMB 
maps as well. For this purpose, besides the WMAP raw data, the analysis software is also 
needed for understanding the reasons of these discovered systematics and flnding approaches 
to remove them. It will be helpful if the WMAP data processing package can be also publicly 
opened. 



^The software we used in this work can be found at http://dpc.aire.org.cn/wmap/09072731/release_v2 
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